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ON THE POTTS MODEL PARTITION FUNCTION IN AN 
EXTERNAL FIELD 

LESLIE M. MCDONALD AND IAIN MOFFATT 



Abstract. We study the partition function of Potts model in an external 

(magnetic) field, and its connections with the zero-field Potts model partition 

function. Using a deletion-contraction formulation for the partition function 

f5 Z for this model, we show that it can be expanded in terms of the zero-field 

C^ partition function. We also show that Z can be written as a sum over the 

^~^ spanning trees, and the spanning forests, of a graph G. Our results extend to 

VO Z the well-known spanning tree expansion for the zero-field partition function 

f\] that arises though its connections with the Tutte polynomial. 

o 

s,J 1. Introduction 



Let G be a lattice or, more generally, a graph. The g-state Potts model partition 



with Hamiltonian 



Ci function at temperature T in an external field H is given by 

(Te5(G) 

> 

Cn (1) hia)^-J J2 5{<^^.<^J)-H J2 '5(1' '^*)- 

^ {v,,v,}eE{G) v,<£V(G} 

C^ Here the Vi are the vertices of G and {vi,Vj} its edges; the cr^ G {!,..., g} are 

r~^, the classical spin variables at u^; J is the spin-spin coupling; /3 = 1/(kT); and the 

f^ sum is over all states a of G. (A fuller description of notation can be found in 

^~~l Section [2J) The Potts model has not only proved to be invaluable to the study 

^^ of phase transitions and critical phenomena in physics (see, for example, [H 1211 

^ |36]), but has also found applications as widely varied as tumour migration, foam 

behaviours, gas absorption, and social demographics (see, for example, [18, 221 [131 
[211 [2S1 ISI] ) • In addition, it admits significant connections with the mathematical 
C^ areas of graph theory (see, for example, [3 HH [251 [M]) and quantum topology 

(see, for example, [HI [37]). Here we are interested in connections between the 
Potts model and graph theory, particularly with graph polynomials. It is well- 
known that these connections have been very fruitful, with the extra generalization 
provided by them leading to significant advances in both combinatorics and physics 
(for example, [3 H [H [HI HT] ) . 

Connections between the Potts model and graph polynomials have traditionally 
focussed on the zero-field Potts model {i.e., when H = in Equation Q), whose 
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2 L. MCDONALD AND I. MOFFATT 

partition function we denote by Zz^ro- For example, the classical, and well-known, 
relationship between the Potts model the Tutte polynomial, T{G;x,y)^ given by 

(2) Z,,ro{G- q, v) = gfc(G)„l^(G)hfe(G) J.(G; (g + v)/v, V + 1), 

where v = e^'' — 1, requires the absence of an external field. (See [3 HI EHl [32] 
for expositions of this connection.) Although this relationship has resulted in valu- 
able interactions between combinatorics and statistical mechanics, (particularly for 
problems involving computational complexity and the study of the zeros of these 
polynomials O HSl |22l |28l |34| ) , there is a need to extend the theory, and in partic- 
ular the deletion-contraction formulation of Zzero provided by its connection with 
the Tutte polynomial, so as to incorporate the external fields (i.e. when H y^ 0) re- 
quired by applications of the Potts model to physics and to other areas that require 
models with external influences. 

Only very recently have external fields been investigated in the context of graph 
polynomials. In the special case of the Ising model (which is essentially the Potts 
model with q = 2) an interesting new polynomial has been found in |33| . This poly- 
nomial has a deletion-contraction relation and captures approximating functions, 
but does not specialize to the Tutte polynomial, nor does the deletion-contraction 
extend to non-constant magnetic fields. Some weighting strategies have also been 
developed to study non-zero magnetic fields in a graph colouring framework (see 
[ini [m |29] ) , with an accompanying graph polynomial. These papers show that the 
Potts model partition function with an external field term does not have a tradi- 
tional deletion-contraction reduction. It follows that Z{G) can not be written as 
an evaluation of the Tutte polynomial T{G). To overcome this difficulty, in [T^], an 
extension of the concept of the contraction of an edge of a vertex weighted graph 
was introduced. (In terms of the Potts model, the vertex weights are used to record 
contributions of the external field H, as described below.) Using this definition of 
contraction, an extension of the Tutte polynomial, called the V-polynomial was 
introduced, and it was shown that the Potts model partition function Z{G) with a 
Hamiltonian given by (II]) is an evaluation of the V-polynomial. A number of con- 
sequences of this result were described in [12] , including a Fortuin-Kasteleyn type 
representation for Z{G) (see also dUl HH ISSl I3S]) j a deletion-contraction relation, 
and computational complexity results. This work unified an important segment of 
Potts model theory bringing previously successful combinatorial machinery, includ- 
ing complexity results, to bear on a wider range of statistical mechanics models. 

Here we exploit the combinatorial machinery provided by this new formulation of 
Z(G) as a graph polynomial V with a deletion-contraction-type relation to obtain 
a number of new expansions for Z{G). Below, we will in fact use a slightly more 
general (but still standard) Hamiltonian than that exhibited in Equation (II]). (We 
favour the more familiar (flj) in the introduction for expositional purposes.) For 
our results, we allow variable spin-spin coupling and site-dependent external fields. 
This extra generality is to allow for further applications of our results, such as to 
the random field Ising model. Below we use the Hamiltonian from Equation (|4|, 
and we will denote the corresponding partition function by Z^xtiG), noting that it 
includes Z{G) as a specialization. 

We expand the Potts model partition function in an external field, Zext{G), in 



terms of the zero-field Potts partition function, ZzeroiG) (see Theorem 3.1). It is 
well-known that the Tutte polynomial, T(G]x,y), can be written as a sum over 
spanning trees. This spanning tree expansion of the Tutte polynomial gives rise 
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to a spanning tree expansion of ZzeroiG). We extend this result, in Theorem 3.2 



to a spanning forest expansion and a spanning tree expansion for the Potts model 



partition function Z{G) external field. Furthermore, we show, in Theorem 3.3 
that Zf^xt (G) admits the spanning tree expansion. This spanning tree expansion for 
Zext (G) extends the spanning tree expansion for the zero- field Potts model partition 
function Zzero{G) that arises through its connection with the Tutte polynomial 
shown in Equation (pi). We conclude by specializing our results to several other 
common types of Hamiltonian from the physics literature. In addition to these 
results, a secondary purpose of this paper is to introduce the new tool of the V- 
polynomial to the physics literature. 

This paper is structured as follows: our notation is described in Section [2J Full 
statements of our main results on Z^xtiG) are given in Sectional with their proofs 
deferred until Sect ions Island [6] The main tool used here is the realization of Z^xtiG) 
as an evaluation of the V-polynomial, which is described in Section [4j 

This work is based on L.M.'s Masters Thesis which was supervised by LM. 

2. Notation and terminology 

2.1. The g-state Potts model. We begin by recalling some basic information 
about the Potts model. Let G be a graph and consider a set {1,2, ... ,q} of q 
elements, called spins. A state of a graph G is an assignment of a single spin to 
each vertex of the graph. Thus, if the vertex set of G is V{G) = {vi, . . . , w„}, then 
a state of G is a function a : V{G) — >■ {1, . . . ,q}. We let S{G) denote the set of 
states of G. 

The interaction energy may be thought of simply as weights on the edges of the 
graph. In physics applications, the interaction energies are typically real numbers. 
However, in our context of graph polynomials, they may simply be taken to be 
independent commuting variables. We will denote the interaction energy on an 
edge e = {vi,Vj} by Je = J^j = Jv,,vy 

The zero field Hamiltonian is 



(3) h{a) = ~ 2_^ J,j(5(cr,,crj), 

{h3}eE{G) 

where ct is a state of a graph G, where ai is the spin at vertex Vi, where S is the 
Kronecker delta function, and each edge {i,j} is assigned an interaction energy Jij. 
Here we will use a common generalization (see, for example, [TJ [71 [5^1 [IHl 155]). 
of the Hamiltonian shown in Equation ([l]). 

Definition 2.1. Let G be a graph. Assign to each edge e = {i, j} an interaction en- 
ergy J^ = Jij", and assign to each vertex Wi a /leW vector Afi :— {Mi^i,Mi^2, ■ ■ ■ ,M^i,q)- 
Then the Hamiltonian of the Potts model with variable edge interaction energy and 
variable field is 

9 

(4) ^(ct) = - Yl J^,A<^^.<^3)-' Yl Y^^i'^^^'^'^''^- 

{i,j}eE Vi£V(G)a=l 

Note that in this model the external field contributes Mi^^i G C to the Hamil- 
tonian for each vertex Vi with spin a^. The field vectors just record the possible 
contributions M^j- at each vertex. (The formal tool of field vectors is used here 
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Figure 1. Deletion and contraction of an edge in a vertex- 
weighted graph. 



in order to naturally describe the deletion-contraction relation for the Potts model 
partition function in an external field.) 

Note that the Hamiltonian of Equation (|3| may be recovered from the general 
Hamiltonian of Equation Q by taking all the field vectors to be zero. The Hamil- 
tonian also subsumes other frequently studied forms of the Hamiltonian from the 
literature, including (fTl). We refer the reader to [H] for details. 

Regardless of the choice of Hamiltonian, the Potts model partition function is 
the sum over all possible states of an exponential function of the Hamiltonian: 

Definition 2.2. Given a set of q spins and a Hamiltonian /i, the q-state Potts 
model partition function of a graph G is 

Z(G)- J2 e-^'-^'^ 
(Te5(G) 

Furthermore, we let 

• Zzero{G) denote the Potts model partition function with the Hamiltonian 
of Equation (Is]); 

• ZextiG) denote the Potts model partition function with the Hamiltonian of 
Equation Q. 

2.2. Vertex-weighted graphs. We use standard notation for graphs. E{G) and 
V{G) denote the edge set and vertex set, respectively, of a graph G. If A C 
E(G), then n{A)^ ^(^); ^ind k{A) are, respectively, the nullity, rank, and number 
of components of the spanning subgraph^ {V{G),A), of G with edge set A. 

Here, a vertex-weighted graph consists of a graph G, with vertex set V{G) = 
{i'i,i;2, . ■ . ,fn} equipped with a weight function w mapping V{G) into a commu- 
tative semigroup. The weight of the vertex Vi is the value w(fi). We will often 
write UJi rather than uj(vi), and x for an indexed set of variables. An edge-weighted 
graph G is a graph equipped with a mapping 7 from its edge set E{G) to a set 
7 :— {7e}ee-E(G)- We use the standard convention 7 : e i-7> 7e, for e e E{G). In 
addition, at times we will write "7^" instead of "7ei"- In this paper, we assume 
that all graphs are edge-weighted and vertex-weighted. 

If G is a vertex- weighted graph with weight function w, and e is an edge of G, 
then G — e is the graph obtained from G by deleting the edge e and leaving the 
weight function unchanged. If e is any non-loop edge of G, then G/e is the graph 
obtained from G by contracting the edge e and changing the vertex weight function 
as follows: if Vi and Vj are the vertices incident to e, and v is the vertex of G/e 
created by the contraction, then w(u) = i^{vi) -\-uj{vj). See Figurell] Loops are not 
contracted. 
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3. Expansions for the Potts partition function 

3.1. Relating Zea;t and Z^^^^ via connected partitions. Let G be a vertex- 
weighted graph with vertex set V{G), and let vr — (Vi,...Vfe) be a partition of 
V{G). Each vertex set Vi induces a subgraph Gi of G. (The subgraph Gi consists 
of the vertex set Vi and all of the edges of G for which both ends lie in y^.) We 
say that a partition tt of G is a connected 'partition if each induced subgraph Gi, 
for i = 1, . . . , fc, is connected. We let V{G) denote the set of connected partitions 
of V{G). We let C(7r) denote the set of subgraphs {Gi,...,Gfe} induced by a 
connected partition tt. 



fe 



Finally, we let [x*] denote the coefficient operator, so that if P = X^i^i 



OjX- 



then [cc'jP = Oi. 

We can now state the first of our main results. This result expresses the Potts 
model partition function with site dependent external fields, in terms of the zero 
field Potts model partition function. 

Theorem 3.1. Let G be a graph equipped with a field vector Mi = (Mi.i , . . . , Mi^q) £ 
C at each vertex Vi. Then the Potts model partition function admits the following 
expansion: 

7reP(G) HeC(TT) 

where, for each connected partition tt, C^tt) is the set of subgraphs of G induced by 
tt; X{tt) = Xv^ ■ ■ ■ Xv^ for Xm — X]q=i e^^^" , and Vi is the sum of the weights, 
Mi , of all of the vertices in Vi . 

We defer the proof of this theorem until Section |6] 

3.2. Spanning tree and spanning forest expansions for Z^xt- Our other two 
expansions for Z^xt require some additional notation on trees, forests and edge 
activity. 

A spanning forest of G is an acyclic spanning subgraph. We let J^{G) denote 
the set of spanning forests of G. If G is a connected graph then a spanning tree of 
G is a connected spanning forest. Following a standard abuse of notation from the 
theory of the Tutte polynomial, if G is a graph with components Gi, . . . , G^g)i 
then we say that T is a spanning tree of G if the restriction of T to each component 
of G is a spanning tree of that component, i.e. TO Gi is a spanning tree of Gi for 
each i = 1, . . . , k{G). We let T{G) denote the set of spanning trees of G. 

Suppose T e T{G) is a spanning tree of G. If e € T, then the cut defined by e 
is the set {/ e -B(G)|(r- e) U / S Tifi)}; and if e ^ T the cycle defined by e is the 
unique cycle in T U e. If, in addition, the edges of G have been given an ordering 
{ei, . . . , Cm}, then edge e G T is internally active if it is the smallest edge in the cut 
defined by e, and is internally inactive otherwise; and an edge e ^ T is externally 
active if it is the smallest edge in the cycle defined by e, and is externally inactive 
otherwise. Activities are defined analogously when F G .7-"(G) is a spanning forrest 
of G (simply replace ''T{G)" by ''F{Gy\ and "tree" by "forest" in the above). 
We will let I^(G, T), Xi{G, T), £a{G, T) and Si{G, T) denote the set of internally 
active, internally inactive, externally active, and externally inactive edges of G with 
respect to the spanning tree T or spanning forest T . When the graph G is clear 
from context, we will sometimes write Xx{T) instead of Xx(G,T), and do similarly 
for the other sets above. Also note that since 2x(G, T) and Ia{G, T) are both sets 
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of edges of the spanning subgraph T, we can regard Tx{G, T) and IaIG, T) as sets 
of edges of T itself. This allows us to write expressions such as T/2i(G, T). 

Our second main result is a spanning forest expansion for the Potts model par- 
tition function with site dependent external fields. 

Theorem 3.2. Let G be a graph equipped with a field vector Mi ~ {Mii, . . . , Mj „) € 
C? at each vertex Vi. Then the Potts model partition function with Hamiltonian (HI) 
admits the spanning forest expansion 










where X{F) = Xp^ ■■■Xp^,, for Xm — X]q=i^ ^^"^ '^"■'^ ^i '■^ ^^^ ■''^"^ ''/ ^^^ 
weights, Mi, of all of the vertices in the i-th component, Fi, of the spanning forest. 

The proof of this theorem appears in Section l6] 

Our final main result is a spanning tree expansion for the Potts model partition 
function with site dependent external fields. It is well-known that the Tutte poly- 
nomial can be defined as a sum over spanning trees (see, for example, the exposition 
in [5j). Using this spanning tree expansion for the Tutte polynomial together with 
Equation ([2]), one can obtain a spanning tree expansion for zero- field Potts model 
partition function with constant edge interactions: 



(5) Z,,ro{G;q,eJ" -1) = 



^fc(G)(-g/3j_-^yy(G)|-fc(G) y- | g + '^ ' ~M (g/sjy^-^^^)! 



T<^T(G) 

(Note that this spanning tree expansion can be extended to a spanning tree expan- 
sion for the zero-field Potts model partition function with variable edge interactions 
using Equation ([7| and Corollary 5.2 below.) Here we extend the expansion shown 
in Equation (jsl) to Zext{G): 

Theorem 3.3. Let G be a graph equipped with a field vector Mi = {Mii, . . . , Miq) £ 
at each vertex Vi . Then the Potts model partition function admits the spanning tree 
expansion 
(6) 

ZeAG)^ J2 f n (e^''=-l)|f n <^"A ZeAT/Ix{G,T)). 

T(iT(G) \eeIx(G,T) / \eG£^(G,T) / 

The proof of this theorem appears in Section [5] 

4. The ^-polynomial and the Potts model 

4.1. The Tutte polynomial and the Potts model. There are seminal relations 
between the Tutte polynomial, its extensions, and the zero-field Potts model par- 
tition function Zzero{G). The multivariate Tutte polynomial, Zt{G), (see |28j and 
[30' ) of a graph G with edge-weights 7 := {7e}ee-E(G) is defined by 

Zt(G;0,7):= E ^'^^^ H ^- 

A<ZE(G) eeA 
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It is well known that the Potts model partition function with the zero-field Hamil- 
tonian and variable edge-interaction energies of Equation ([3]) can be recovered from 
the multivariate Tutte polynomial: 

(7) Z,,ro{G) = Zt[G- q, {e^'^ - l}eeB(G))- 

(See |13j for the early stages of this theory and later exposition in, for example, 
jH [51 ISl 1311 [31].) Furthermore, if the edge-interaction energies are constant with 
Je — J, for each edge e, then the Potts model partition function can be recovered 
from the Tutte polynomial as in Equation (12]). As mentioned in the introduction, 
this connection has resulted in a remarkable interactions between the areas of com- 
binatorics and statistical mechanics, particularly for computational complexity and 
the study of the zeros of these polynomials. (See, for example, [3 [THl HHIIHI [33! for 
surveys of these results.) The equivalence of the Tutte polynomial and the Potts 
model partition function assumes the absence of an external field, and Z^^t can 
not be expressed in terms of the (multivariate) Tutte polynomial. (Some weight- 
ing strategies have been developed to study non-zero fields in a graph colouring 
framework. See [IHlinilSni.) 

4.2. The V-polynomial and the Potts model. In [12], EUis-Monaghan and 
Moffatt used a definition of contraction that incorporates vertex weights to assimi- 
late a Hamiltonian of the generic form of Equation Q into the theory of the Tutte- 
Potts connection. To do this, a new graph polynomial, called the V-polynomial, 
was introduced. The V-polynomial generalizes Noble and Welsh's VF-polynomial of 
[191 , and extends the Tutte polynomial by incorporating vertex weights and adapt- 
ing contraction to accommodate them. It was shown that the variable field Potts 
model partition function (with its many specializations) is an evaluation of the V- 
polynomial, and hence is a polynomial with a deletion-contraction reduction and a 
Fortuin-Kasteleyn type representation. 

Definition 4.1 ([E]). Let 5 be a commutative semigroup, let G be a graph 
equipped with vertex weights uj := {uji} C S and edge weights 7 := {"fe}, and 
let X — {xk}keS be a set of commuting variables. Then the V -polynomial, V(G) — 
V{G,uj;x,j) e 'Z[{le}eeE{G)A^k}kes], of the vertex- and edge-weighted graph G, 
is defined recursively by: 

(1) V(G) = V(G - e) + j^V{G/e), if e is a non-loop edge of G; 

(2) V{G) = (7e -I- l)y(G - e), if e is a loop; 

(3) V{E,n) — YliLi x^- , liEm consists of to isolated vertices of weights wi, . . . , w„ 

We will use the following state sum expansion for the V-polynomial later. 
Theorem 4.2 ([12]). V(G) can he represented as a sum over spanning subgraphs: 

ACE(G) eeA 

where Ci is the sum of the weights of all of the vertices in the i-th connected com- 
ponent of the spanning subgraph (V{G),A). 

The main result of [12 was the recovery of Zext{G) as an evaluation of the 
Tutte-type graph polynomial V{G), extending the classical connection between 
T, Zt and Zzero, and providing a deletion-contraction reduction for the partition 
function Zf,xt{G). 
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Theorem 4.3 ( 12J). Let G be a graph equipped with afield vector Mi = (Mi.i, . . . , Mi^q) g 
C? at each vertex Vi, and let h{a) he the Hamiltonian given in Equation (Hi). Then 

Zext{G) = V {G,uj; {XmJmgCi, {e'^ '•' - l}{ij}e£'(G)) , 
where the vertex weights are given by a;(wi) — Mi and, for any M — (Mi, . . . , Mq) € 

We note that Theorem |4.3| unifies an important segment of Potts model theory 
and brings previously successful combinatorial machinery, including complexity re- 
sults, to bear on a wider range of statistical mechanics models. Furthermore, The- 
orem 



4.3 gives a Fortuin-Kasteleyn-type representation for Zext{G) and shows that 



the partition function Zf,xt{G) is in fact a polynomial. See |12j for details. 
5. Spanning trees and the Potts model 



In this section we prove Theorem |3.3[ which gives a spanning tree expansion for 
the Potts model partition function Zext{G). This result generalizes the spanning 
tree expansions for Z^ero that come from the Tutte polynomial T{G) and the mul- 
tivariate Tutte polynomial Zt{G). The expansion for the partition function Ze^;* 
will follow from the spanning tree expansion for the V-polynomial given in the 
following theorem. 

Theorem 5.1. Let G be a graph equipped with vertex weights u) := {^i} and edge 
weights 7 ;= {je}- Then 

V{G,u;x,j)^ J2 I n 7e) I n (7e + l)| V(T/Tx(G,T)). 

TeTiG) \e6lx(G,T) / \ee£^(G,T) / 

Proof. For ease of reading, we set 

a(G,r):=( n 7e)( n (7e + l)). 

eeIx{G,T) ee£A{G,T) 

We will prove the theorem by induction on the number of edges of G. 

The assertion is easily verified when G has no edges. 

Now suppose G has ordered edges {ei, . . . , e„i}, rn > I, and assume the assertion 
holds for all graphs with fewer than m edges. Note that the order of the edges of 
G induces an order of the edges of G — e^ and of G/e^- We will use these induced 
orders throughout the proof. We will consider three cases: when Cm is ordinary, a 
bridge, or a loop. 
Case 1: If Cm is ordinary, then 

(8) ^ a{G,T)ViT/IiiG,T)) 

TeT(G) 

= J2 a{G,T)V{T/Ix{G,T))+ ^ a{G,T)V{T/Ii{G,T)). 

TeTiG) T<£T{G) 

e^l^T e„eT 

Focussing on the first term on the right-hand side of pi, we have that T e T{G) 
with em ^T ii and only if T e T{G — em). Since e™ is not a loop and is the largest 
edge. Cm G £i{G, T). In addition, for 1 < i < m, we have that Ci € £^(G, T) if and 
only if Ci G £j\{G — em,T), (this is since the cycles defined by e^ in G and G — e^, 
with respect to T, are identical). We also have that e^ € Xx{G,T) if and only if 
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Ci e Ii{G — em,T), (this is since the ordinary edge e™ is the largest edge and so 
is never the smallest edge in the cut defined by e^ in G and in G — e„i)- When 
e™ i T, we then have £:^(G, T) = £^(G - e„, T), and Ii{G, T) = Ii{G - e^,T). 
Thusa(G,r) = a(G-e„,r) and r/Ii(G,r) ^ T/2i(G - e™,T). 

Focussing on the second term on the right-hand side of (Is]), we have that 
T £ T{G) with em G T if and only if T/cm € T{G/em)- Also, since e™ is not 
a bridge and is the largest edge, e™ G Ii{G,T). In addition, for 1 < i < m, we 
have that e^ G £^{G,T) if and only if e^ G £^{G/ern,T/em), (this is since the 
cycles defined by e^ in G and G — e,„, both must contain an edge smaller than e„i). 
Furthermore, we have that e,; G Ti{G,T) if and only if e^ G Xx(G j CnnT j e^) ^ 
(this is since the cuts defined by e^ in G and G/cm, with respect to the two 
trees, are identical). When e^ G T, we then have £x{G,T) = £x(G/em,T/em), 
and Ix{G,T) = Ix{G / e^.T / e„^) U {e™}. Thus q;(G, T) = 7,„a(G/e„j,r/e„) and 
r/Zx(G,T) = r/(2i(G/e,„,r/e™) U {e™}) = (T/em)/Zi(G/e„,r/e„). 

Using the above observations, we can then write (Is]) as 

Y, a{G~e„,,T')V{r/Ix{G~e„,,T')) 

T'eT{G-e^) 

+ 7m Y. «(G/e™,T")F(r7li(G/e„,r")), 

T"er(G/e„) 

which, by the inductive hypothesis, is equal to V{G— em) +'ymV{G/em) = V{G), 
as required. 

Case 2: If e^ is bridge, then V'(G) — V{G — €„) + 7my(G/e„j), which, by the 
inductive hypothesis, is equal to 

Y a{G-em,T')V{T'/Ix{G~em,r)) 

T'(£T{G-e^) 

+ 7m Y a{G/em,T")V{T"/Ix{G-em,T"))- 

T"eT{G/e„,) 

Since e„i is a bridge, we have that T is a spanning tree of G if and only if T — e^ 
is a spanning tree of G — e™. Also, T is a spanning tree of G if and only if T/cm 
is a spanning tree of G/em- 

Also since Cm is a bridge, it is the unique edge in its cut, and so e,„ is not 
in £j^(G,T) nor in Xx{G,T). Furthermore, since Cm is the largest edge, checking 
the activity of an edge e^, 1 < i < m, in G with respect to G is the same as 
checking the activity of e^ in G — e™ with respect to G — Cm, which is also the 
same as checking the activity of e^ in G/e„i with respect to G / e„i- Thus we have 
£ji,{G,T) = £j^{G-em,T-em) = £^^(G/e,„,T/e„i) andIi(G,T) =Tx{G-em,T- 
em) =Tx{G/era,T/em)- Thus a(G,T) = a(G-e,„,T-e™) = a(G/e,„,r/e„i) and 
r/Zi(G, T) - (T - e„0/2i(G - e„., T - e„0 = (r/e„0/2i(G/em, T/e„0. 

Using this, we may write the above sum as 

Y ["(G - ^m, T -em)V {{T - e„j)/2i(G - e,„, T - e^)) 

TeT{G) 

+ 7™ a{G/e„i, T/e„i) V ((r/e„i)/2i(G/e„, T/e™))] 

= ^ a(G,r)(V([T/Zx(G,r)]-e„0 +7„F([T/Ix(G,r)]/em)), 
T£r{G) 
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which, via deletion-contraction, is just X^TsTfG) '-'^(^'^) ^ (-^/-^2:(G, T)), as re- 
quired. 

Case 3: Finally, if e™ is a loop, then V(G) = {jm + 1)V{G — e„,), which, by the 
inductive hypothesis, equals 

(9) (7m + 1) Y. «(G-e„,T')F((G-e„0/2x(G-e™,T')). 

T'<£r{G~e,^) 

Clearly, T G T{G) if and only if T e T{G — em)- Since Cm is a loop, it is externally 
active with respect to every spanning tree T of G. Also, for 1 < i < m, the 
activity of e^ in G with respect to T is identical to the activity of e^ in G — Cm 
with respect to T. Thus a(G, T) = (7™ + l)a(G - e^, T - e^) and T/Ii(G, T) = 
TjXxiG - em, r - e™). 

Using this, the expression (lo]) can be written as 

^ a(G,r)y(r/Zi(G,r)), 

Ter(G) 

and the result follows. 

D 

It was shown in [12 that the multivariate Tutte polynomial Z{G) can be recov- 
ered from the V-polynomial by setting Xi =^ 9 for each i: 

(10) F(G,w;x,;=0,7)-^t(G;0,7). 



We can use this fact together with Theorem |5.1| to recover Traldi's spanning tree 
expansion for the multivariate Tutte polynomial from [301 (see also [S]). We will 
use this expansion in the proof of Theorem |6.2| 



Corollary 5.2 (Traldi '30'). Let G be a graph equipped with edge weights 7 :~ {7e}- 
Then 

ZT{G;en)^o'^^^ J2 U ^e li iie + e) II (7e + i), 

TeTiG) eGXz(G,T) eeXMG,T) eG£^(G,T) 

where k{G) denotes the number of components of G. 
Proof. By Equation p^ and Theorem |5.1l 



Zt(G;0,7)= E n ^- n (7e + l)^T(T/Xi(G,T);0,7). 

Ter(G)eeIz(G,T) ee£^(G.T) 

However, T/Ii{G,T) is just a forest with fc(G) components whose edge set is in 
bijection with Ia{G,T), so Zt {T/Ii{G,T)) = 9''^^^ neei^(G,T)(^ + le), and the 
result follows. D 



We note that the well-known spanning tree expansion for the Tutte polynomial 

^|I^(T)|„|£^(T)| 
T£T[G) 



T{G;x,y)= ^ a:l---"y 



can easily be recovered from Corollary |5.2[ In addition, Theorem |5.1| gives spanning 
tree expansions for Noble and Welsh's polynomials U{G) and W{G) from [TU]. 
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Remark 5.3. Observe that the proof of CoroUary |5 . 2| indicates why the expression 
V {T/Xx{G,T)) appears in the spanning tree expansion for the V-polynomial, but 
Zt does not appear in the expansion for the muhivariate Tutte polynomial (and 
similarly for the Tutte polynomial): while the values of the Tutte polynomial and 
the multivariate Tutte polynomial on trees are easily computable, this is not the 
case for the V-polynomial. 



Proof of Theorem \3.3\ By Theorem |4.3[ we have 

Z(G) = V {G,uj; {XMJMeC^, {e^^= - l}eeB(G)) 
The result then follows by applying Theorem [Sj 



6. Spanning forests and the Potts model 



In this section we prove the expansions stated in Theorems |3.1| and |3.2[ Our 
approach to both of these theorems is to find a corresponding expansion for the 



V-polynomial and then to apply Theorem 4.3 Both of the expansions for the V- 
polynomial and the proof of these expansions are extensions of results due to Noble 
and Welsh on the VF-polynomial in |19j . Our first result shows that the polynomial 
V(G) admits an expansion in terms of the multivariate Tutte polynomial Zx- 

Theorem 6.1. Let G be a graph equipped with vertex weights a; :— {uji} and edge 
weights 7 := {7e}- Then 

v{G,oj-xn)^ J2 ^W n [l']{ZT{H;qn)), 

7rG-P(G) HeC{7r) 

where a;(7r) := Xc-^ ■ ■ ■ Xc^, and Ci is the sum of the vertex weights in the i-th block 
of the partition tt. 



Proof. We begin with the state sum expansion of V{G) from Theorem 4.2 

y(G;a3,7)= ^ Xc,Xc^ ■ ■ ■ Xc^^^^Y[-/e- 

ACE(G) e<£A 

Letting n(^) denote the partition of V{G) whose blocks Vi consist of the vertices 
in the connected components of the spanning subgraph with edge set A, we have 

(11) V{G;X,-f)^ Y^ Xc,Xc2---Xc,^^^Y[-fe= Yl ^('^^ Yl II'^^' 

ACE{G) e<EA Tre'P(G) A:n(A)=7r eSA 

where the sum is over all connected partitions tt of V{G). 

Suppose that the partition tt e ^^{G) has blocks Vi, ■ ■ ■ ,Vk and the connected 
components of the subgraph induced by tt are Gi, . . . ,Gk where Vi = V{Gi). Then 

E n^e- E n > 

A:n(A)=7r eSA AiU---uAfcC_B(Gi)u---UB(Gfc) eeAiU---uAfc 

(VijAi) connected 
AiCE{Gi) 

k k 

= n E n> =l[iq']Z{G.;q,l)= n {'l']Z{H;qn)- 

i=l Ai<ZE(Gi) e^Ai i=l HeC{7T) 

(Vi,Ai) connected 



Substituting this into Equation (11) gives the result. D 
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Proof of Theorem \3.1\ The result follows immediately from Theorems |4.3| and |6.1[ 
and Equation ([t]). D 



Using the sum over connected partitions from Theorem |6.H we will obtain an 
expression for V{G) as a sum over spanning forests of G. We will use this spanning 
forest expansion to prove Theorem |3.2[ 



Theorem 6.2. Let G be a graph equipped with vertex weights a; := {uJi} and edge 
weights 7 := {7e}- Then 

viG,io;x,j)= Y. x{F) n > n (i+>) 

FeJ^iG) eeE{F) ee£A(G.F) 

where x{F) := a;^ ---Xc^, where ciis the sum of the vertex weights in the i-th 
component of F. 

Proof. By Theorem |6.1[ we have 

(12) F(G;a;,7)= ^ x{n) [] [g^] (Z^li/; 9,7)) • 

TTev(G) Hec(TT) 

Also, by Corollary |5.2[ when H is connected, 



= E n > n 7e n (i+>)= e n ^e n (i+»- 

TeT{H) eeXx{H,T) eeX^(H,T) e££j^{H,T) TeT[H) eeE{T) ee£A{H,T) 

Substituting this into Equation ( |12[ ) gives 

V{G;xn)= J2 ^W n E n > E i^ + le) 

TreV{G) HeC{Tr)Ter{H) e<£E{T) ee£^(T) 

= E ^w E n > n (i+>) 

n(F)=7r 

= E ^(^) n > n (!+>)' 

FeJ^iG) eeE(F) ee£A{F) 

as required. D 

We can now prove the second of our main results. 



Proof of Theorem 3.2 The theorem follows immediately from Theorems 4.3 and 6.2 



D 

7. Applications to common models 

In this final section we provide some specializations to a few other common 
models. 

The first example is used in models in which a system with an external field in 
which one particular spin (the first, without loss of generality) is preferred (see for 
example the surveys [51 136) V In the standard applications to magnetism in models 
of magnetism in an external field /3 = 1/{kT), where T is the temperature of the 



ON THE POTTS MODEL PARTITION FUNCTION IN AN EXTERNAL FIELD 13 

system, and where k = 1.38 x 10^^'^ Joules/Kelvin is the Boltzmann constant. Here 
again, variable interaction energies and field values are allowed. 

Theorem 7.1. Let G be a graph in which a complex value Zi is associated to each 
vertex Vi, let h be the Hamiltonian 

(13) h{a) = - ^ Ji^jS{ai,aj) - ^ ZiS{l,ai), 

{»j}e-B v,ev(G) 

and let Z denote the Potts model partition function arising from this Hamiltonian. 

(1) //, for zeC, X,=eP^ +q-l, then 

Z(G)= Y. ^W n [q']{Z.ero{H;q,{e^-'^ -lUEiH))) . 

TxeV{G) HeC{n) 

where, for each connected partition -k , X{'k) = X^-^ ■ ■ ■ X^^, where Cj is the 
sum of the weights, Zi, of the vertices in the j-th block of the partition n, 
andC{-K) is the set of subgraphs of G induced by tt. 

(2) //, for zeC, X^^eP"-- +q-l, then 



z{G)= y: x(F)('n(e^"-i))( n ^ 



ee£^(G,F) 



where X{F) = X^ • ■ ■ X^^, where Cj is the sum, of the weights, Zi, of the 
vertices in the j-th component of the spanning forrest F. 
(3) Z admits the spanning tree expansion 

^(G)= E f n (e^'''-i)i( n ^^'^\z{Tiix{G,T)). 

TeTiG) \eGlx(G,T) 



Proof. The theorem follows easily from Theorems |3.1[ |3.2[ and |3.3 




upon noting 



that the Hamiltonian (13 1 can be recovered from the Hamiltonian (4|) by setting 
Mi = {zi, 0, . . . , 0) for each i, and that then each vector Mi is described by z^ and 
so we may reindex the expressions using complex numbers. D 



The following corollary contains the restriction of Theorem 7.1 to models with 
constant interaction energies and a constant, non-zero field. 

Corollary 7.2. Let G be a graph, h be the Hamiltonian 

(14) hia)^-J Y 5{a,,a,)~H ^ 5{l,a,), 

{i,j}eE v,£V{G) 

and let Z denote the Potts model partition function arising from this Hamiltonian. 
(1) We have 

Z{G)^ Y X{7t) n [q']{Z.ero{K;q,eP-'~l)), 
TTeV(G) Kec{TT) 

where, for each connected partition it, X{'k) — Ilili (*^^ +9^1); where 
kin) is the number of blocks in n, and \Tri\ is the size of its i-th block. 
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(2) We have 

z[G)= y: x{F)(i[ic^''^-i)]( n ^''^ 

FeF(G) VeGF / \ee£A(G,F) 

where, X{F) = 0^=1 v'^^'^^^'''' +9^1); where Fi is the i-th component of 
the spanning forrest F. 

(3) Z admits the spanning tree expansion 

^(G)= E f n («''■'= -1))( n ef'A Z{T/MG,T)). 

TeT{G) \eeXx{G,T) J \eGf^(G,T) / 

Our final specialization is to random field Ising model (RFIM) which is used to 
study disordered states. It has a random field in that the z^'s are randomly chosen 
local magnetic fields that each affect only a single site. In [T^], it was shown that the 
RFIM (which is essentially the q = 2 Potts model) can be expressed as an evaluation 
of the V-polynomial. The RFIM is also particularly interesting from the point of 
view of graph polynomials as there have been two very recent studies of the Ising 
model from a graph theoretical perspective. [2] treats the Ising model with constant 
interaction energies and constant (but non-zero) magnetic field as a graph invariant 
and explores graph theoretical properties encoded by it. |33) on the other hand, 
creates a new graph polynomial, which has a deletion-contraction reduction for 
non-loop edges, and which is, up to a prc-factor and change of variables, equivalent 
to the RFIM. 

The RFIM takes spin values in { — 1, 1}. We will let r denote a state for the 
RFIM, which is a map t : V{G) -^ { — 1, +1}. As usual we set r^ := T{vi). Also we 
will let T{G) be the set of states for the Ising model. 

Theorem 7.3. Let G be a graph with a vertex weight z^ G C associated to each 
vertex Vi . The RFIM given by 

h{T) = -J Y. ^»^J ~ J2 ^»^^' ""'^ ^('^) = J2 e-'3''(^) 

{i,j}<£E i<£V{G} r£T(G) 

admits the following expansions. 
(1) We have 

Z(G) = c-^''^^) Y a;(F) (e^^^ - 1)1^(^)1 \{ e^^' , 

FeJ^{G) ee£_^(G,F) 

where for any z ^ C, Xz = c^^ + c''^,- ri{G) = J\E{G)\ + iJ2i(£V(G) ^ijO.'^d 
x{F) := Xc-i ■ ■ ■ Xck, for Cj the sum of the vertex weights in the j-th compo- 
nent of F. 



(2) Wehave, forri{G) = J\E{G)\+3J2 



Zi 



ieV{G) ^«' 



Z(G) = e"^''^'^^ > (e^P-^ - l)l^^(-^^l g^A'>^|t.A('J')l+^';C-^'Ai('J'))_2'(T/Ii(T)). 



Proof By Theorem 6.3 of [ig, Z{G) = e'^''^^^'>V {G,uj; {x^jzeC, e^^'^ ~ l), where, 
for any z & €., Xz = e^^ -|- e^^, and a;(wi) = z^. The first item follows by applying 
Theorem |6.1[ and the second item follows by applying Theorem |5.1[ 

D 
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